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Integrable Many-Body Systems via Inosemtsev Limit 



The Inozemtsev limit (IL) or the scaling limit is known to be a procedure applied to the elliptic 
Calogero Model. It is a combination of the trigonometric limit, infinite shifts of particles coordinates 
and rescalings of the coupling constants. As a result, one obtains an exponential type of interaction. 
In the recent paper it is shown that the IL applied to the sl(N,C) elliptic Euler-Calogero Model 
and the elliptic Gaudin Model produces new Toda-like systems of N interacting particles endowed 
with additional degrees of freedom corresponding to a coadjoint orbit in sl{n, C). The limits cor- 
responding to the complete degeneration of the orbital degrees provide only ordinary periodic and 
non periodic Toda systems. We introduce a classification of the systems appearing in the sl(3, C) 
case via IL. The classification is represented on two-dimensional space of parameters describing the 
infinite shifts of the coordinates. This space is subdivided into symmetric domains. The mixture 
of the Toda and the trigonometric Calogero-Sutherland potentials emerges on the low dimensional 
domain walls of this picture. Due to obvious symmetries this classification can be generalized to 
the arbitrary number of particles. We also apply IL to sl(2, C) elliptic Gaudin Model with two 
marked points on the elliptic curve and discuss main features of its possible limits. The limits of 
Lax matrices are also considered. 



1 Introduction. 

1. Inozemtsev limit Let us consider a system of two interacting particles with the 
hamiltonian: 



where v and — v their momenta in the center of mass frame, u is the difference between their 
coordinates, m is the coupling constant and E%{u, r) is the Eisenstein function defined on 
the complex torus T 2 with moduli r. It is related to the Weierstrass functions as follows: 



Inozemtsev limit (see [3], [4], [8]) is a combination of the trigonometric limit Im(r) — > oo 
with the infinite shifts of coordinates and a scaling of the coupling constant. This procedure 
can transform the above hamiltonian into the following one written in the shifted coordinates 
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E 2 {z,t) = P {z,t) + 2C,{\). 



U = u- 



i/m(r) and with the rescaled coupling constant m: 



H = v 2 + m 2 {e u + e~ u ) 



(2) 
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In a similar manner the IL may be applied to the system of N interacting particles: 

N N 

H(vi,Ui) = +m 2 ^E 2 {u l -Uj) (3) 

i=l i^j 

The limit leads to the periodic Toda system (non periodic case is also possible): 

N N 

ffK £/*) = ]>>? + 2m 2 eUi ~ Ui W 

i=l i,j=i+l;j=l,i=N 

Thus IL describes the connection between the Elliptic Calogero-Moser and the Toda models. 

2. Purpose of the paper The purpose of the paper is to apply the IL to different 
kinds of generalizations of the elliptic Calogero model and to find out what sort of systems 
appear in the limit depending on its parameters. 

We are going to apply the IL to the following models: 

1. The sl(N, C) elliptic Euler-Calogero Moser model (ECM): 

L k i(z,z) = v k S k i - pkkSklEi(z) - p k i(l - 5 k i)e 27Ti ^ Uk '<P(u k i,z) (5) 
where p^ e sl(N, C) and Poisson brackets are: 

{vi,Uj} = Sij, {pij,pu} = PkjSu - puSkj, {pij,Vk} = 0, {pij,Uk} = (6) 

N 

H( Ui , Vi) = ^ v\ + PijPji E 2 ( U iJ ) ( 7 ) 
»=1 ijtj 

The integrable system is defined on the symplectic factor sl(N,C)/ /{diag SL(N,C)}. In 
this case it is given by p k k = and some gauge fixing condition (see section 2). In the 
following the degrees of freedom related to pij on the reduced space are called the orbital. 
The hamiltonians Hk,i are defined by the Lax matrix Lij(z 7 z) G sl(N, C), z e T 2 : 

N 

Tr(L k (z,z))=Y, E ^ H ^ 

1=0 

where Ei(z) <~ are the elliptic functions (see Appendix A). 

2. The elliptic Gaudin model (EG) (see [6]): 



No 

{L G )kl = V k Skl - y^ J {Pa)kkSklEl(z - X a )- 
a=l 

-(I -Ski) E e 2wilZ ~ Xal - r * a)ukl (Pahi*(u k i,z - x a ) 
where {p a )ki S sl(N,C). The Poisson brackets are: 

{Vi,Uj} = Sij, {(Pa)ij, (Pb)kl} = Sab{{Pa)k]5il - (Po)«^fej)) 
{{Pa)ij,U k } = 0, {(Pa)ij,V k } = 0, a = l...N 



(8) 



(9) 



3. Notations For making the limits it is convenient to consider the torus T 2 : C/ (2wiZ + 
2o) 2 Z), w\ = —in, Im{w2) = 0, r = 

The trigonometric limit is produced by UJ2 — > oo. Bar is used for the rescaled variables 
and function \ f° r the powers of the exponents: m = me XmU! , \m = x( m )- 

The shifts of coordinates arc introduced as follows: 



Uij — Uij + ,fijUJ 2 



(10) 



4. Results In the sl(3, C) ECM model the quadratic hamiltonian follows from (6): 

3 

h 2 ,o = J2 v i~ 2aE ^ u ^) - 2 bE 2 (u 23 ) - 2cE 2 {u 31 ) 

where a = p\ 2 p 2 i, b = P23P32, c = P31P13 are the orbital variables on the reduced phase 
space (the last one is d — pi2P23P3i)- 

Let us write down the possible quadratic hamiltonians, obtained from the above one via 
the IL, modulo permutations of the particles and degenerations (some terms can vanish in 
the limit): 

3 

1. h 2 ,o = v l - 2 ae- Ul2 - 2be~ U23 - 2ce~ U31 (11) 

i=\ 
3 

2. h 2 . = J2^- 2de Ul2 - 2be~ U23 - 2ce~ U31 (12) 

i=l 

3 

3. h 2 ,o = ^ v i - 2a(e _c/l2 + e Ul2 ) - 26e" c/23 - 2ce~ U31 (13) 

4. h 2 , = J>? - \d-^ ~ 26e"^3 _ 2 - ce u 3 , (14) 
i=i 2 

5. h 2 ,o = Y,vf- \d— ^ - 2b(e- U23 + e U23 ) - 2c(e t *» + e"^) (15) 

i=l 2 

' f-f 2 sinh 2 ^ 2 sinh 2 %i 2 sinh 2 V 7 

2=1 2 2 2 

The first and the last one are a direct generalizations of the periodic Toda and the Suthcrland- 
Calogero systems correspondingly. All others content an unusual dependence on the coordi- 
nates. 

The classification can be given by the following figure (fij are parameters of the shifts of 
coordinates: tiy = Uij + u 2 fij and thus fij + fjk + fki = 0): 

X. //23 = 

/12 = 2- 



/23 = 1 



/l2 = 




/31=0 /31 = -1 /31 = -2 &gA 



On fig .1 different kinds of the systems correspond to symplexes of different dimensions 
and situations (modulo permutations of a,b,c): {1} and {III} - two types of two dimensional 
symplexes, which content the limits to hamiltonians of the first (11) and the second (12) 
types (domains {II} and {IV} are equivalent to in sense of the permutations); lines 



{fij = 2Z} and midlines {fy — 1 + 2Z} content intervals which correspond to the mixture of 
the Sutherland and the Toda (14) and cosh(J7jj) terms (13); all vertices - zero dimensional 
symplexes correspond to fifth (15) and sixth (16) types of the hamiltonians. 

The considering of the cubic hamiltonians provides more complicated structure on the 
fig.l. by adding a dual lattice (see fig. 2 in Section 4). 

The resultant picture (unification of fig.l and fig. 2) appear to describe the classification 
of possible limits of brackets between a, b, c and d. 

Fig. 1,2 have obvious symmetries and thus the classification can be generalized to an 
arbitrary number of particles. In the case of N interacting particles all possible limits of 
hamiltonians and Poisson brackets correspond to different types of domains in the N — 1 
dimensional hyperplanc of variables fij . 

There are more possible limits in the EG because the diagonal part of the Lax matrix 
(9) content dependence on {p a )u even on the reduced phase space and there are additional 
parameters (x a —Xb) which can be also shifted in the limit. In this case resultant hamiltonians 
content a direct dependence on the parameters of the limits (see Section 7). 



2 EG and sl(N,C) ECM. 

The Elliptic Gaudin Model is an example of the Hitchin system [2] (see also [5] and [7]) and 
was introduced in [6]. The Lax matrix is given by (8). 
We can restrict algebra (9) on the submanifold: 

^2(p a )kk=0, k = l...N ( 17 ) 

a 

since it corresponds to the fixing of the moment map fi p by the adjoint action of the diagonal 
SL(N, C) subgroup on p k i and the fact pkk are Poisson commute with the Tr(L(z) n ). By 
fixing the adjoint action one reduces the number of independent variables in (p a )ki by TV — 1. 
Note that the restriction (17) is in agreement with the simple requirement on holomorphic 
function: 

res{(L G (z)) kk } = ^ 

simple poles 

The hamiltonians appear as the coefficients in the equation describing the spectral curve: 

det(L(z) + m) = (19) 

Another way is to decompose Tr(L m ) on the basis of functions {E\(z — x a ), E 2 (z — 
x a ), E' 2 (z — x a ), ...}. Then the coefficients in the decomposition are again hamiltonians (see 
[6],[5]). 

The sl(N, C) ECM corresponds to N = 1. In this case (17) is written as follows: 

{» P )i=Pii=0 i = l...N (20) 

and thus the Lax matrix is: 



(L) w = v k 6 k i - e 2m *-> kl p k i<f>(u kl ,z) 

k^l 



(21) 



Using appendix B we obtain to the following hamiltonians in this case: 

<L(zf >=H 2t0 + H 2 . 2 E 2 (z), (22) 
< L(z) 3 >= #3,0 + H 3t2 E 2 (z) + H 3t3 E' 2 (z), (23) 



N 



#2,0= -PijPji E 2{Uij)), (24) 

AT 

#2,2 = 51 Pyft'' ( 25 ) 



JV 

#3,0 = K 3 - 'iviPijPjiE 2 {u ij ) - PijPjkPki[E2{uij)E-i{ujk) + E 2 (u j k)E 1 (uij)+ 

i,3,k=l (26) 

+^E 1 (u k i)(E 2 (u i j) + E 2 (uj k )) + ^E^Uij - Ujk){E 2 (uij) - E 2 (u jk ))}), 

N 

#3,2 = ^ (3 v iPijPji + PijPjkPki[E\(Uij) + E^Ujk) + Ei(u ki )]), (27) 

1 w 

#3,3 = X] PioPjkPki- (28) 
»,j=l 

Note that i?2,2 and #3,3 are the Casimir functions of the algebra (6). 

It should be also mentioned that the exponential factors in the Lax matrices (5), (8) 
are not necessary for computing the hamiltonians since tiy + Uji = 0, Uij + uj k + u k i = 0. 
However these factors plays an important role in the limits of the Lax matrix itself (see 
Section 8). 



3 s/(3,C) ECM. 

Let us compute the dimension of the phase space of the sl(N,C) ECM. In the Pij space 
a number of the independent variables equals (N — 1)(N — 2) since the integrable system 
is defined on the orbit {pij / / Diag SL(N,C)} with the values of the Casimir functions 
fixed, where DiagSL{N 1 C) is the Cartan subgroup in the SL(N, C). In the space (u, v) £ 
T*Bunc,T. there are 2(N — 1) such variables. Thus the total dimension of the phase space 
equals N(N — 1) and consequently dim(sl(3, C) ECM) — 6 that is in the agreement with a 
fact there are three integrals of motion in this case: #2,0, #3,0 7 #3,2- 

The orbit {pij/ /Diag SL(N, C)} is given by conditions pa — and some fixation of the 
coadjoint action by the Cartan subgroup. However a certain fixation of the coadjoin action 
(cither a fixation of values of the Casimir functions) may be not useful. Thus we will only 
imply some choice of the gauge fixing condition without any exact one when consider not 
invariant quantities. 

Poisson algebra (6) between nondiagonal variables in TV = 3, N = 1 case looks like: 

{Pl2,P23> = -Pl3, {P32,P2l} = ~P31, 

{P31,f>12} = -P32, {P21,f>13} = -P23, (29) 
{P23,f>3l} = -P21, {P13,P32} = ~Pl2- 

We may do not care about fixing the action of Diag SL(N, C) if choose variables which are 
invariant with respect to it: 



a = pi 2 p 2 i, b = p 23 p 32 , c = pi3p 3 i, d = p 12 p 23 p 3 i 



(30) 



Note that a, b, c and d are the combinations of pki sufficient for introducing the hamiltonians 
(24 — 28). From (29) we obtain the brackets between the invariant variables: 



{a, b} = {b, c} = {c, a} = d - 
{a, d} = ab — ac, {b, d} = be — ab, {c, d} — ac — cb 



(31) 



Note also that there is no singularity in ^ since it equal to P21P32P13- Hamiltonians are 
given by (24) - (28): 



h 2 ,o = v\ + vl + v\ - 2aE 2 (u 12 ) - 2bE 2 (u 23 ) - 2cE 2 (u 31 ) 



'2.2 



2(a + b + c) 



h 3 ,o = vf+v^+ vl + 3v 3 aE 2 (u 12 ) + "ivxbE^u^) + 3v 2 cE 2 (u 31 )+ 

3 abc E' 2 {u 12 )E 2 {u 23 ) - E' 2 (u 23 )E 2 {u 12 ) 
2 [ d ' E 2 (u 12 ) - E 2 (u 23 ) 

h 3 , 2 = 3(t>i + v 2 )a + 3(vi + v 3 )c + 3(v 2 + v 3 )b+ 

+ 3(rf - ^)[^l(«12) + £i(«23) + El(«3l)] 

3 . , abc. 



(32) 
(33) 

(34) 
(35) 



(36) 



Here h 2 . 2 and ft.3,3 are the quadratic and the cubic Casimir functions of the algebra (31). 
The spectral curve is written as follows: 



det(L(z) + /1) = 
/1 3 - lfi{E 2 (z)h 2t2 + h 2 . ) + \{E' 2 {z)h 3t3 + E 2 (z)h 3 . 2 + h 3fi ) = 



(37) 



4 Degenerations of elliptic functions. 

1. Formulas used for computing IL in quadratic hamiltonians 

Let us put x = x n + <toj 2 , where x n is a new (shifted) coordinate, 
\a\ < 2, < (3 < 1, u) 2 -> 00 
Below we right down the main non vanishing order in the limit uo 2 — > 00. 

l. a = 0: ^coth(^) 



1 ^ x 

= 9 X, «>th(- -fcwa) 



fc— — OO 



2. < cr < 1 : ^ +e -*n-™ 2 

3. cr = 1 : - - 2e" W2 sinh(a;„) 

4. cr > 1 : - - e Xn-(2-*)u2 

2 



k—oo 



e x+ 2 kui 2 -\- g- 



- 2 ku> 2 — 2 



-CTW2 + 2o) 2 2 



(38) 



(39) 



1. (3 = 0, 0<a<2 
1.1 a = 0, (3 = 0: -- 



4 sinh 2 ^ 



(40) 



1.2 (3 = 0, < a < 1 : e "^-x n 

1.3 /3 = 0, cr = 1 : e - u '(e- Xn +e Xn ) 

1.4 /? = 0, 1< cr < 2 : e -(2-^2+x„ 
2. /3 > 0, < a < 2, 

2.1. cr = /? < 2 - cr : e"*" 

2.2. 2-cr = (3<cr: e x " 
2.3 cr = /3 = 2-cr: e" x " + e x ~ 

From (40) one can simply obtain the quadratic hamiltonians (11) — (16) and (4) given in the 
Introduction. For example (4) appears if one put Uj = Uj + 2jui 2 -^ m = rheTi^ 2 (see [8]). 

2. Formulas used for computing IL in cubic hamiltonians 

1. < cr < 1 : - e -™»-*« 
2. o = 1 : e-"*(e x " -e" x ") 



1 y COSKI - k. 2 ) 



(41) 



3. 1< cr < 2 : e -(2-fV2+x n 



Let us consider the limits of following expression: 



A _ 1 E 2 (u 12 )E 2 (u 23 ) - -^2(^23)^2(^12) 
2 £*(ui 2 ) -£ 2 (u 23 ) 



(42) 



It is contented in (13^ given by (34). To analyze this expression one should take into account 
terms of the second non vanishing order in the approximate formulas for E 2 (z) and E' 2 (z). 
The result is represented on fig.2: 



1 E 2 {u 12 )E 2 (u 23 ) - E' 2 (u 23 )E 2 {u 12 ) 

2 E 2 (u 12 ) - E 2 (u 23 ) 



A 

A_ : -e 



e Ul2 e ~ (2-/12 )W2 



-t^l2 g — /l2^2 



B_ : -e 



-c/23. 



)w 2 
/23^2 



(43) 



_ e -C/31g-(2+/3l)w2 



Areas ^4/,, B'_, B' + and C(_ are shifted on the hole period with respect to A_, A + , 
B + and C_ correspondingly. 

There arc two types of boundaries on fig. 2. The first corresponds to the case when the 
boundary belongs to one of the axis (f\ 2 € 2Z or f 23 G 2Z or /31 e 2Z). All other lines are 
contented in the second type. They belong to the dual lattice. In fig. 2 one can see a hexagon 
A' + B'_C' + A'_B' + C- . The sides of this hexagon content the boundaries of the first type and 
the lines coming through the centers of the hexagons content the boundaries of the second 
type. Let us write down an example of limits of A corresponding to the first type of the 
boundaries: 

c+ n {/ 23 

A r. < * (44) 

A+n{f 2 3 = 0}:e u ^e- 



0} : -e-^e-^l-coth^p) 



-(2-/ 12 ) ( l + C oth^) 



The boundaries of the first type correspond to a combination of the trigonometric limit and 
the IL and thus content the mixture of the Toda and the Calogero-Suthcrland potentials. 

As for the boundaries of the second type the values of A from the neighboring domains 
are of the same order and thus the potential is equal to the sum of those from the neighboring 
domains. (All the same is valid for the zero dimensional boundary.) 



3. sl(N,C) case In the N particle case one should consider tr(L 2 ),...,tr(L N ). The most 
complicated potential appears as a free term in the expression: 



$(iti2, z)$>(u 2 3, z)...$(u N1 ,z) 



(45) 



It is possible to generalize fig. 2 on higher dimensions. In sl(4, C) case the space of the 
parameters is three dimensional (f\ 2 + f 23 + ^34 + /41 = 0). The whole space is subdivided 
into regular tetrahedrons with a three dimensional structure inside that correspond to the 
fig. 2 on the two dimensional walls. The case of the arbitrary N can be obtained by analogy. 



5 Limits of brackets for s/(3,C) case. 

We are going to find non-singular limits of the brackets (31) and all possible limits of the 
hamiltonians (32 — 36) for each of them. 

The idea is to multiply if it is necessary a certain hamiltonian by the exponent with the 
power canceling the highest power in the hamiltonian. 



(46) 



This procedure saves the involution property of the hamiltonians but certainly does not 
guaranty the independence of hij. The sufficient condition for the independence of the 
hamiltonians is a saving of the kinetic terms ^ vf and ^ vf in h 2t o and Z^o correspondingly. 

Another (equivalent to the previous) way of studding the limits comes from considering 
the spectral curve (37). 

H 3 - ^(E 2 (z)h 2 , 2 + h 2 . ) + ^{E 2 (z)h 3 , 3 + E 2 (z)h 3 . 2 + h 3fi ) = 

One may require to have a non singular spectral curve in the limit. We will sec it is true for 
all limits if x(^3,o) = x(h 2 ,o) = 0. 

One can also make the rescaling of the spectral parameter: 

/i = fie x » (47) 

However in this case the kinetic terms vanish in 1i 2 _q and h 3 ,q and this may lead to the lost 
of independence of the hamiltonians. To save the independence one may resale momenta 

Vi = Vie x " (48) 

But this can be done only if there is a factor e~ Xfi in the canonical form: 

3 
i=l 

In this sense the rescaling of the momenta does not provide interesting consequences. We 
consider the case x M = 0, xiMfi) = 0, x(h 3 . ) = 0. 

The purpose of the section is to find possible (non-singular) limits of the brackets (31): 

{a, b} — {b, c} — {c, a} — d — , {a, d} = ab — ac, {b, d} = be — ab, {c, d} — ac — cb. 

and thus find out whether there exist limits with non-trivial dynamics of rescaled a, b, c, d: 

a = e x ^ 2 a, b = e XbUJ2 b, c = e x ^ 2 c, d = e XdUJ2 d. (50) 

The non singularity conditions for the rescaling of the brackets (31) are: 

Xa < Xd, Xb < Xd, Xc < Xd, 

(51) 

Xd<Xa+Xc, Xd<Xa+Xb, Xd<Xb + Xc- 

It is important to notice that it is impossible to satisfy the above non-equalities if some 
X<0. 

Non triviality of the rescaled brackets means an exact equality besides any of the inequal- 
ities (51). Let us write down all such possibilities modulo permutations of a, b, c: 



Xa < Xd, Xb < Xd, Xc < Xd, Xa+Xb = Xd, Xa+Xc> Xd, Xb + Xc > Xd 

1 : { (52) 
{a, 6} = d, {b, c} = 0, {c, d} = 0, {a, d} = 0, {6, d} = 0, {c, d} = 

Xa < Xd, Xb < Xd, Xc < Xd, Xa+Xb = Xd, Xa+Xc = Xd, Xb + Xc > Xd 

2 : { (53) 
{a, 6} = d, {b, c} = d, {c, a} = 0, {a, d} = 0, {b, d} = 0, {c, d} = 

Xa < Xd, Xb < Xd, Xc < Xd, Xa ~\~ Xb — Xd, Xa + Xc Xd, Xb ~\~ Xc — Xd 

3 = { _ _ _ _ _ _ (54) 

{a, b} = d } {b,c} = d 7 {c,a} — d 7 {a, d} — 0, {b, d} — 0, {c, d} = 



4: ^ 



5 : < 



Xa = Xd, Xb < Xd, Xc < Xd, Xa+Xb> Xd, Xa+Xc> Xd, Xb + Xc > Xd 

{a, b} — 0, {b, c} — — ^S-, {c, a} = {a, d} — 0, {b, d} — — ab, {c,d} = ac 

Xa = Xd, Xb = Xd, Xc < Xd, Xa+Xb> Xd, Xa+Xc> Xd, Xb + Xc > Xd 

dbc 



{a,b} = 0, {b,c} = {c,a} = -- 



d ' 



6 : ^ 



7 : < 



{a, d} — ab, {b, d} = —ab, {c, d} = ac — cb 

Xa = Xd, Xd = Xd, Xc = Xd, Xa+Xb> Xd Xa + Xc > Xd, Xb + Xc > Xd 

{a, b} = {b, c} = {c, a} = 

{a, d} — db — dc, {b, d} = bc — db, {c, d} = ac — cb 

dbc 
T' 



Xa = Xd, Xb = 0, Xc = Xd, {a, b} = {b,c} = d- 



{c, a} — 0, {a, d} — {d, c} — —ac, {b, d} — be — ab 



: < 



0<Xa<Xd, 0<Xb< Xd, Xc = Xd, Xa + Xb = Xd, {a, 6} = d - 



(55) 



(56) 



(57) 



(58) 



(59) 



{b, c} = {c, a} = 0, {a,d} = —ac, {b,d} = bc, {c, d} = 
Let us note that the spectral curve unites h^ : 2 and the cubic Casimir /i 3i3 in the limit 



uo 2 — > oo if 



^ 1- 



6 Classification of Limits in s/(3, C)ECM 

The classification consists of possible limits of the hamiltonians (32 — 36) and the limits of 
the brackets (31). 

A notable fact about classification of the limits of the brackets takes place if one require 
non triviality of all three terms in the quadratic hamiltonians in the limit. In this case \a, 
Xb and Xc satisfy some equality coming from /i 2 + /23 + / 3 i = 0. In the domain {/} on 
fig.l it looks like Xa + Xb + Xc = 2 and in {77/}: Xa = Xb + Xc- These equations can be 
either compatible or not with the relations determining some given limit of the brackets. By 
comparing them with the list of the limits of the brackets given in Section 5 we arrive to 
the fact all types of brackets (52) — (59) are distributed on fig. 3 with the midlines drawn on 
fig.l. The domains {77}, {777} and {IV} correspond to brackets 8 - formula (59); midlines 
- type 8, (60) with \d = 1; the sides of the fundamental triangle -type 7, (59); the dual 
lattice figured on fig. 2 split the domain {7} on fig.l by three parts: they correspond to the 
type 1 and 4, (52), (55); and one dimensional boundaries between them - type 2 and 5, see 
(53), (56) and the center of the fundamental triangle - type 3 and 6, see (54), (57). 



Examples Let us write down three examples corresponding to the cases 1, 7, 8 ((52), (58), (59) 
correspondingly) for the limits of the brackets. 



Example 1 

{a, b} = d, {b, c} = 0, {c, a} = 0, {a, d} = 0, {b, d} = 0, {c, d} = 

, 1, 1 , 7 1 1 1 7 

Ji2 — g, J23 — J3i — _ Y2' * a _ 3' * b ~ 4' ^ — 2' ^ d — 12 

h 2 ,o = v{ + w| + vl - 2de- Ul2 - 2be' U23 h 2 . 2 = 2c (60) 

h 3 ,o =v\+vl+vl + 3v 3 de- Ul2 + 3wifee _c/23 - 3de U31 

i ; — n - n dbc - 3- 

fi3,2 + /i3,3 = -3cw 2 - 3— p h 3 . 3 = --d 

Example 2 

{a, b} = {b, c} = d — —*j^> {c, a} = 0, {a, J} = — ac, {6, J} = be — ab, {c, J} = ac 

/l2 = 1, /23 = 0, / 3 1 = -1, Xa = 1, Xb= 0, Xc = 1, Xd = 1 

lg L 

2 sinh 2 (^) 

3,1 , ( 61 ) 



h 2 ,a = vj+vl+ vl - 4acosh(£/ 12 ) - -b — 2 „ - - 4ccosh([7 31 ), /i 2>2 = 2 (a + c) 



(62) 



^3,0 = vl + vl + v 3 + 6v 3 acosh(U 12 ) + ^v x b : " + 6w 2 ccosh(J7 31 )- 

2 sinn 

-3(d - ^){sinh(;7 12 ) + cosh(;7 12 ) coth(^)} 
a ^ 

t „/ \ „, x 3 . — abc, , ,C/ 2 3. 7- 3.- abc, 

h 3 . 2 = 3(«i + u 2 )a + 3(«i + u 3 )c + - (d - ) coth(— ), h 3y3 = -^ d+ ^ 

Example 3 

{a, b} — d~ {b,c} = {c,d} = 0, {d,d} = —ac, {b,d} = bc, {c, d} = 

f 1 f 1 / 1 1 1 i 

Ji2 — 2' '' 23 = 2' = _ *° _ 2' * b _ 2' Xc ~ * d ~ 

h 2 ,o = v{ + v\ + vl - 2de- Ul2 - 2be- U23 - 4ccosh([/ 31 ) 
^2,2 = 2c 

h 3 ,o = vl + vl + vl + 3v 3 ae~ Ul2 + Sv^e" 1123 + 6v 2 ccosh(£/ 3 i) - 3(d- ^) e U31 

t x 3, T abc. T 3,7 abc, 

h 3 , 2 = 3(wi + v 3 )c+-(d- -=-), hs, 3 = ~2^ d+ ~d^ 

7 Limits in si (2, C) EG with two marked points on ellip- 
tic curve. 

Let us consider the following Lax matrix: 

Lm = v k 8 kl + (p 1 ) kk E 1 (z - xi) + (p 2 ) kk E 1 (z - x 2 ) + 

„ . (2 — a;-, ) — (z — aTi ) n . (2 — x 2 ) — (2 — x~ 2 ) (^^) 

+ {Pi)kie 2 ™ " sw $(uh, z - Xi) + (p 2 ) kl c 2 ™ ±=^-""${u,a,z - x 2 ) 



Vl = 


— «2 


« (Pl)ll = 


-(Pi) 


. (P2)ll = 


~(P2) 


f (Pl)ll = 


-(P2) 


1 (Pl)22 = 


~(P2) 



(64) 



22 



Due to (18) we have 

I (Pi)n = -fe)n 

(65) 



Let us use the following notations v = vi, u = u\ 2 and pn = (pi)n Instead of the fixing of 
the action by diagonal SL(2, C) matrix we again choose invariant coordinates: 



(^1)12(^1)21 = ai (^2)12(^2)21 = a,2 (^1)12(^2)21 = a-3 (Pi)n = 04 (66) 
The brackets are easily computed: 

{01,02} = {01,^4} = {a2,«4} = 

{01, a 3 } = {a 2 , a 3 } = 2a 3 a 4 (67) 
{a 3 ,a 4 } = a 3 

To compute the hamiltonians produced by the tr(L 2 ) we need the following identity: 

$(w, «!)$(-«;, z 2 ) = -*(«;, Zl - - #1(22) + Si (to) - £7i(io + «i - z 2 )] (68) 

The hamiltonians are: 

2^2,0 = v 2 - aiE 2 (u) - a 2 E 2 (u)+ 

+a 3 e 2 ^~ X ^ 2 ~ X2) u ^(u,x 2 - x 1 )[E 1 (u + x 2 - Xl ) - £i(u)]+ (69) 

+ gi^ e 2^ ( "- a >:<.' i -* i) « $(U)a;i _ + Xl _ X2 ) _ El (u)], 

03 

7^2,1,1 = 2a 4 t; - 2alE 1 (x 1 - x 2 ) - a 3 e 27ri ~ EL ^ ?2 ^'"^(u, x 2 - Zi)+ 

&1& 2 9-t- ( a 2- 3 =2)-(^l- 3 =l) 

+ _L^ e 2« T _ f -2:2), (70) 

1 7 _ ,2 1 1, _ 1 2 

2^2, 2,1 — Ol + a 4, 2 2' 2 ' 2 ~ a 2 + a 4- 

The notable features of the EG case are: 

1. There is non trivial non diagonal part in the matrix pa that makes hamiltonians and 
algebra more complicated. 

2. There is a parameter x\ — x 2 which can be also shifted. 

3. The resultant hamiltonians directly depend on the parameters of the limit 



8 Limits of Lax matrix. 

Using formulas from the Appendix C one can obtain the approximate expression for the Lax 
matrix (6), -/V = 1, pa = 0: 

L kl = v k S kl - e^^^pki^uki, z) (71) 
It degenerates into the following one: 



Lij « Vi + 

+e -^-fij £ ^-^™2p ije "2Cii±.( co th^- + coth|)(l - e "»+*-**>) » 

« Vi + (72) 
fij > : e-^-f't^^pije-^-™ 2 *" 2 ^ (1 ~ e Ui ' +fi ' UJ2+Zn+a ' JJ2 - 2uJ2 ) 

+ : 



The simplest way to obtain the above formula (besides of the straightforward evaluations 
given in the Appendix C) is to use series representation: 



nez 



(73) 



Let us consider possible limits. 



1. For fij > the requirement of non singularity of when w 2 — > °o is the following: 



-fyv + Cij <o 



fij 



<? + Cij +fij+<r-2<0 



Three exists three non trivial limits for the non diagonal part of : 



w 4 + < 



1 : e~ — 

Ui 



Vij 



2 : e--* a -t» :a ¥ 1L Pi i {-e u *i +x ») 

k 3:e-^ CT -^ £ ^p 4 ,(l- e c/ -+ z ") 
The cases 1, 2 and 3 correspond the following conditions 



2 a 

1 = /ij = -Cij, fij < 2 - cr, Cij < cr(l - 77)- 



2: /i; 



Ci; 



1 - - 

1 2 



fij > 2 - a, Cij < <t(1 - tt)- 



3 : = 2 - cr, Cij = o"(l - -). 



(74) 



(75) 



(76) 



2. For < we have the following non singularity condition for Li 

-fyv + Cij -° < 



-fyv + tij+fij <0 



^i + < 



4 : e 2 CT 4 p y ( e z ~ ) 



5 : e -^ i ^-fe £l 4 £2! -p y .(-e c/l 



{ 6:e-* L <'-*' !n ¥ iL p ij (e-*»-e u «) 
The cases 4, 5 and 6 correspond to the following conditions: 

4 : fij = -Cij - 2, fij < -cr, Cij < 0"(1 - 

cr Z 



5 : 



: /ij = -Cij 1 _ g ' A? > ~ CT > Cij < cr(l - 7^' 



6 : /ij = -cr, Cij = cr(l - -). 



(77) 



(78) 



(79) 



3. Example With the help of the obtained formulas it is possible to obtain the Lax 
matrices for the examples given in Section 6. For the example 3 (see (62)) we have 



L(z) = 



hui2 



pise" 5 ^(e^i _ e -») \ 



P2ie 2 



it/12 



«2 



-P 23 e 2^3 



"3 



(80) 



9 Conclusion. 



In the recent paper it is shown that the IL applied to the sl(N,C) elliptic Euler-Calogero 
Model and the elliptic Gaudin Model produces new Toda-like systems of N interacting 
particles endowed with additional degrees of freedom corresponding to a coadjoint orbit 
in sl(n,C). The limits corresponding to the complete degeneration of the orbital degrees 
provide only ordinary periodic and non periodic Toda systems. We introduce a classification 
of the systems appearing from the sl(3, C) case. The classification is represented on two- 
dimensional space of parameters describing the infinite shifts of the coordinates. This space is 
subdivided into symmetric domains. The mixture of the Toda and the Sutherland-Calogero 
potentials emerges on the low dimensional domain walls of this picture. Due to obvious 
symmetries this classification can be generalized to the arbitrary number of particles. We 
also apply the Inozcmtsev limit to si (2, C) elliptic Gaudin Model with two marked points on 
the elliptic curve and discuss main features of its possible limits. The limits of Lax matrices 
are also considered. 

We are grateful to A.Chervov, A.Marshakov, A.Gorsky, A.Mironov, A.Morozov and 
A.Zabrodin for useful remarks and to A. Levin for many discussions. We especially grateful 
to M.A.Olshanetsky for initiating the work and useful discussions. The work of both authors 
was partially supported by RFBR grant N00-02-16530 and the program for support of the 
scientific schools 00-15-96557. 

10 Appendix A: Basic Formulas and Definitions 

The rest of the formulas given in appendices are borrowed from [5], [8] and [10]. 
The basic element is the theta function: 

oo 

0(*,r) = qie-i{e^ z - e~ iwz ) q n )(l - q n e 2niz )(l - q n e~^ iz ), q = e 2mT (81) 

n=l 

The series representation for the Eisenstein functions is convenient to make the limit u 2 — > oo: 

E 1 (z,t) = d z logd{z,T), £7i(*,r)w- + ... (82) 

z 

E 2 (z,T) = -d z E 1 (z,r), E 2 (z,t) ^ ^ + ... (83) 
Relations to the Weierstrass functions: 

{{z,T)=E 1 (z,T) + 2m(T)z (84) 



p(z,r) = E 2 {z,t) - 2th(t) 



(85) 



where 



In the Appendices B and C we use notations which are convenient for making the limits (see 
[8]): 



r = 



U! 2 
U)l ' 



(86) 



u)\ = — Z7r Im(ui 2 ) = 



11 Appendix B: formulas which are used for computing 
the hamiltonians and their limits 



Let us represent the Eisenstein functions in the form convenient to make the limit ui 2 ^ oo: 

(87) 



k— — oo 

oo 

= - A E 



1 



- fc= -oo Sinh2 (f- fcW2 ) 

1 ~ cosh(| - ku 2 ) 
4,^ sinh 3 (f-fcw 2 ) 
The following important expression appears in the Lax matrix: 

+ z)tf'(O) 



It has a pole at z — and 



The well known identities are: 



res| z=0 $(u, z) = 1 



?;)$(— it, u) = £ 2 (u) - -E 2 (u) 



(88) 
(89) 

(90) 



(91) 



& u {u,v)=${u,v)(E 1 {u + v)-E 1 {u)) 



And consequently 



Ei(u + v) + Ei(u — v) — 2E 1 (u) 



E 2 {u) -E 2 (v) 

The following identity is necessary to find the cubic hamiltonians: 
/ \E 2 {u)E' 2 {u) \ 

2-d(v — w)"&(w — u)"&(u — v)-d(u + v + w) 



det 



1 E 2 (v) E' 2 (v) 
V 1 E 2 (w) E' 2 (w) J 



It can be rewritten in the following way: 

*(«, zMv , zm-u - v , z) = -\e> {z) - ^w g[;j:g[g - 

1 £ 2 (i;)£ 2 (ix)-£ 2 (ii)£ 2 (i;) 
2 E 2 (u)-E 2 {v) 

Using formula (93) we have 



-£ 2 (i/) 
E' 2 {v)E 2 {u) - E' 2 (u)E 2 (v) 



2E 1 {u + v)- 2E 1 (u) - 2E 1 (v) 

= -Ei(u + v){E 2 {u) + E 2 (v)) + 



E 2 {u) - E 2 {v) 

+E 1 (u - v)(E 2 (u) - E 2 {v)) + 2{E 2 {u)E 1 {v) + E^E^u)) 
Thus formula (95) can be written in the following form: 

$(u, «)*(-u- v,z) = -Ie^(z)+E 2 (z)[E 1 (u) + E 1 (v)+E 1 (-u-v)]- 

-E 2 {u)Ei{v) - E 2 (v)E 1 (u)+ 

+ ^E 1 (u + v)(E 2 (u) + E 2 {v)) - ^E x {u - v){E 2 (u) - E 2 (v)) 



(92) 
(93) 



(94) 



(95) 



(96) 
(97) 



(98) 



12 Appendix C: formulas which are used for computing 
the limits in the Lax matrix 



tf( ZjT ) = q h-^{e^ z - e-" z ) f[(l - q n ){\ - q n e 2mz )(l - q n e- 2 

71=1 



q = e 2 " T (99) 



where £(z) is the Riemann zeta function 

t?(u,t) wi _i 



i?'(0,t) 7T 



3 12 JT s ^ n ( 7r { u ~ nT )) 



(100) 



(101) 



n— — oo 



= — e jj£ sinh( — \- nujz) 



tf'(0,r) 

oo 

-2ie^ sinh | JJ (e 2 + e~ z - e 2 ™ 2 - e" 2 ™^) 



(102) 



-2i • sinh - JJ(-1 



a z-2nu 2 I — z — 2nuJ 2 _ — 4raa;2 



2i-sinh(-)(l-e 



z — 2u 2 g — z — 2w 2 \ 



If p~ > and > than 

OLU2 OL02 



{\ _ e u+z — 2L0 2 _ e — u — z — 2uj 2 ^ 



1 - e 



u+z— 2o;2 



(103) 
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